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Abstract 

The radial part of the wave function of an electron in a Coulomb 
potential is the product of a Laguerre polynomial and an exponential 
with the variable scaled by a factor depending on the degree. This note 
presents an elementary proof of the orthogonality of wave functions 
with differing energy levels. It is also shown that this is the only 
other natural orthogonality for Laguerre polynomials. By expanding 
in terms of the usual Laguerre polynomial basis an analogous strange 
orthogonality is obtained for Meixner polynomials. 



1 Introduction 

The time-independent Schrodinger equation for one electron subject to a 
Coulomb potential is 

h 2 c 
-—Ai/f (x) - —V (x) = Ei> (x) , 
2m \\x\\ 

for certain constants m,c (see Robinett 0], Ch. 18). Rewrite the equation 

as 

k 

Alp (x) + Tj— tti/j (x) = — Xip (x) , 
\m\ 



'During the preparation of this paper the author was partially supported by NSF grants 
DMS 9970389 and DMS 0100539. 
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where A; is a positive parameter and A denotes the rescaled energy level; let 

x € M. N (usually N = 3, of course) and ||x|| = (X/^Li 2 ^) • We separate 
variables by considering solutions of the form ip (x) = cp(\\x\\)pi (x) where p\ 
is a harmonic homogeneous polynomial of degree I = 0, 1, 2, . . . (the angular 
momentum quantum number). Then (ft must satisfy the equation 

d 2 . . N + 21- 1 d , \ fk 



dr i^ (') + ~ ^ (r) + I - + A J ^ (r) = 0. 



Using spherical polar coordinates for R we see that (ft must be square- 
integrable for the measure r 2 + dr on {r : r > 0} (we restrict our at- 
tention to the point spectrum). Furthermore solutions for the same value 
of I but different energy levels must be orthogonal for this measure, by the 
self-adjointness of the equation. The topic of this note is a direct elementary 
proof of this orthogonality. It will also be shown that these functions and 
the standard set form the only possible (in a sense to be specified) orthogo- 
nal sets of scaled Laguerre polynomials. Then by expanding the functions in 
terms of the usual Laguerre basis an analogous orthogonality is obtained for 
Meixner polynomials. This orthogonal set incorporates a parameter depen- 
dent on the degree and is a set of eigenfunctions of a self-adjoint second-order 
difference operator resembling equation (Q). 

In equation ([[]) set A = — fi 2 (the bound energy levels are negative) and 
substitute <\> (r) = e~^ r g (2 fir), then g (s) satisfies 



n 



The integrable solutions are obtained by setting ^^1 + ^ — ^ — 2/ 
for some n = 0, 1, 2, 3, . . . then (ft (r) = L" (2fi n r) exp (—fi n r) where fi 
2(a+2n+i) an< ^ ot = N + 21 — 2 ] see (H) below {(ft is not normalized). The 

energy levels are — fi? n (specifically — ( 4 ( n +; +1 ) ) f° r N = 3). 

These radial Laguerre-type functions continue to be useful in the re- 
search literature: Beckers and Debergh JjJ] used them for supersymmetric 
properties of the parastatistical hydrogen atom, Blaive and Cadilhac 0] 
derived an expansion in powers of fi n (for fixed I), Dehesa et al Q found 
approximations for entropies which involve integrals of the logarithms of the 
Laguerre polynomials, Xu and Kong || constructed shift operators for the 
radial Hamiltonian. 
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The Pochhammer symbol is defined by 

n 

(a) n = H(a + i-l),aeR, 
i=l 

and the hypergeometric (2-^1) series is 

2 F 1 (a,b;c;x)=jr " J ' ^ V . 

i=o l ^ J ' 

The basic facts about the Laguerre polynomial of index a and degree n 
that will be needed are: 

L° (x) L« (x) x^e^dx = 5 mn r (a + 1) ( " + J )n , (2) 

X -J~? +(a + l- 1)^ + ^1° (x) = 0, (3) 
ax z ax / 

(n + 1) L£ +1 (x) - (a + 1 + 2n - x) L° (x) + (n + a) L^_ x (x) = 0. (4) 

Their proofs and further information are in Szego ||], Ch.5. The author 
gratefully acknowledges a useful discussion with R. Yahez. 



2 Orthogonality relation 

Fix a parameter a > — 1. The scaling factor k appearing in equation (ffl) is 
set equal to 1. 

Definition 1 For n = 0, 1, 2, . . . and x £ R Zei 
n (x) = L° ( — — — — ■ ) exp 



a + 2n + 1 / V 2(a + 2n + l 



Theorem 2 For m, n = 0, 1, 2, . . . f/ie following orthogonality relation holds: 
/ <£ m (x) (x) x Q+1 dx = ^ (a + 2n + l) a+3 T (a + 1 + n) . 
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The proof is composed of two parts. First, the case m = n follows from 
the orthogonality relation @ and the three-term recurrence (Q) in the form 

xhl (x) = - (n + 1) L£ +1 (x) + (a + 2n + 1) L a n (x) + (a + n) L a n _ x (x) , 

and thus xL® (x) 2 x a e~ x dx = (a + 2n + 1) T (a + 1 + n) /n!. In the in- 
tegral f °° </> n (x) 2 x a+l dx change the variable to s = a+ 2n+i • 

Second, fix m ^ n and in the integral J* °° m (x) (x) x Q+1 dx change 
the variable to 

x ( 1 1 



2 \a + 2m + l a + 2n + 1 

and let 

a + 2n + 1 



a + m + n + 1 



then the desired integral equals ( ^—h — — Li^ijl — !Li — 1 j times 

a+m+n+1 / 



L% l (f3s)L«((2-(3)s)s a+1 e- s ds. (5) 



o 



We will show this equals zero by expressing the two polynomials in terms of 
LJ +1 (s), which simplifies evaluation of the integral. Recall the identities 

K (x) = K +1 (x) - Ig\ (x) , (6) 

L S n (bx) = £ {5+ . i y}^ V (1 - b)^ L] (x) , (7) 
i=o ln 

which hold for arbitrary x, 6 E R, 5 > — 1 and n = 0, 1, 2, . . . . Then (with 
5 = q + 1 in (|)) 

L« (Ps)=L^ 1 ((3s)-L%-_\(Ps) 
m ( a + 9 + 7") 

j'=o v 



x ((1 - 0) (a + 1 + m) - (m - j)) 

-2 + j) m 
(m - j)\ 
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and 

L« ((2 -(3)s)= L a n +1 ((2 -IS) a)- L a n + _\ ((2 - /?) s) 
n („ + 4- i) 

= e fr-^r 1- ' ( 2 - # ^ - i ) n_1 " j L r 1 « 

x ((/3-l)( Q + l + n)-(n-j)) 

= E (a + , 2 1 J V W (2 - # (/? - 1)n_1 " J (j - m/3) L T X (s) • 

Now multiply the two expressions together and integrate with respect to the 
measure s a+1 e~ s ds on {s : s > 0}. By the orthogonality (||) with parameter 
a + 1 integral (g|) equals T (a + 2) times 

(a + 2) m _ 1 (a + 2) ri _ 1 



m! n! 

min(m,n) 



(1 - (3) m - L {fl - I)"" 1 (8) 



(a + 2),j! V(l-/3) 2 y 



i=o 



The formula holds when one of m, n = by the convention (a + 2)_ 1 = ;^pj • 
Let 

_ W-2) =1 1 
7 (I-/?) 2 (I-/?) 2 

then 

(j — n (2 — /?)) (j - mm = j (a + 1 + j) H mn - j (a + m + n + 1)) , 

7-1 

since (3 {2 — (3) = The second line of the sum (|8|) can be written as 

7 2 mn(a + m + n+l) 
(7-1) (a + 2) 
which equals van^^ times 

(7 - 1) iF\ (1 - m, 1 - n; a + 2; 7) + 2-^1 (— m, — n; a + 2; 7) 

7(a + m + n + 1) 

— 2 Fi (1 - m, 1 - n; a + 3; 7) ; 

a + 2 

and this expression is zero by the contiguity relations for the 2-Fi-series (by 
the straightforward calculation of the coefficient of 7^ for each j). This 
completes the proof of the theorem. 



mnj 2F1 (1 — m, 1 — n; a + 2; 7) H 2-F1 (— m, — n; a + 2; 7) 

7-1 
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3 Uniqueness 



We show that the only possible orthogonal sets {L° (7 n x) exp (— ^7 n x) }^1_ 
for the measure x a+K dx on {x : x > 0} are the two cases (1) k = 0, -y n = 70 
for all n; (2) k = 1, 7 n = 7o a + 2n+i • By rescaling x we may assume 70 = 1. 
The choice of the basis functions is motivated by their general utility and 
precise asymptotic formulae. Let (j) n (x) = (j n x ) ex P ( — \lnx) (with each 
7„ > 0) and define the inner product 

(f,g)=T(a + K+ I)" 1 / f(x)g (x) x a+ *dx, 

Jo 

for suitable functions /, g. By using the same method as in Section 2 (for- 
mula (H)) the inner product (<f) m , 4> n ) is expressed as a multiple of 

/•oo 

/ L« (/3s) L« ((2 - 0) s) e- s s a+ *ds, 
Jo 

where 8 = 27 ™ . The evaluations are straightforward algebraic calcula- 
tions and were carried out in the computer algebra system Maple™, which 
was also used to find the polynomials q2, <?3 and their resultant. 

3.1 The case k ^ 0,1 

We show that the equations {4>q, (pi) = (0o, 02 ) = ("Pi, <h) = are incompat- 
ible unless k = or 1. The equation (0o, <fii) = implies 71 = a +2n+i • ^ ne 
condition (</>o, <fo) = is equivalent to a certain quadratic equation 52 (72) = 
(not displayed here). The discriminant of equals 16k (a + k + 1) (a + 2). 
Integrability requires a + n + 1 > hence there are real solutions for 72 only 
if k > 0. The conditions (<j>i,(j>2) = and 71 = a +2n+i are equivalent to a 
certain cubic equation q 3 (72) = 0. The resultant of </2)<?3 (a polynomial in 
the coefficients which is zero if and only if there exists a simultaneous root) 
is found to be the product of k (k — 1) by a polynomial in (a + 1) and k 
with all coefficients nonnegative. This shows that q2 (72) = = q3 (72) is 
impossible for n > 0, a > — 1 unless k = or 1. 

3.2 The case k = or 1 

For the case k = the condition {4>o,4> n ) = (for n > 0) and the formula 
(0) with 5 = a show that (/? — l) n = where /3 = 1+ 2 ^ , hence 7n = 1 = 7o 
for all n (the variable 3 is not identical to the one used in Section 2, because 
the scaling parameters 7 n are initially undetermined). 
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For the case k = 1 , (<fo, 01 ) = implies 71 = ; further for n > 2 the 
condition (<pQ,<p n ) = and formulas (^), (0) imply 

09 - if" 1 ((a + 1 + n) (/? - 1) - n) = 0, 

(where /J = ^) with the solutions /? = 1 and fgg^. Thus 7n = 

or 1. But 71 = and the computation of (0i, 0„) with 7 n = 1 yields 

11 1 \ ( U n-i (( .0, (a + 3) (a + 3) n _ 2 

(01, (/>„) = (-1) ra((a + 3)n-l)- — 2 

(a + 1) (a + 2) 

which is nonzero for n > 2. This rules out the possibility that j n = 1 and 
shows that 7 n = q+^+i ^ or an n * s necessary. The sufficiency was proved 
already. 



4 Change of Basis and Meixner Polynomials 

Since there is another obvious orthogonal basis for L 2 (x a+1 <ix), the stan- 
dard Laguerre functions {LjJ," 1 " 1 (cx) e~ cx l 2 : m = 0, 1, 2, . . . }, with some scal- 
ing parameter c, the computation of the transformation coefficients will 
give rise to an infinite orthogonal matrix. The entries will be expressed as 
Meixner polynomials: for parameters 7 > and < c < 1 the Meixner 
polynomial of degree n is given by 

M n (x; 7, c) = 2-F1 (—n, -x; 7; 1 - cT 1 ) , 

with orthogonality relations 

£ { ^M n (x; 7, c) M m (x; 7, c) = , (9) 

£r£ 3:1 (t)„ c n (1 - c) 7 

satisfying the second order difference equation 

c (x + 7) / (x + 1) — (x + (x + 7) c) / (x) + x/ (x — 1) = n (c — 1) / (x) , 

(10) 

and the three-term recurrence 

xM n (x; 7 ,c)= w - ; M n+ i(x; 7 ,c)+ (11) 
c — 1 

J "-^M n x; 7, c + -M n _! x; 7, c 

1 — c c — 1 
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(see the handbook of Koekoek and Swarttouw ||, p. 45). 

Now consider the two orthogonal bases {qo,n ■ n > 0} and {qi t m : m > 0} 
for L 2 ((0,oo) , x a+1 dx) where qo, n (x) = L% +1 ((3x) e"^! 2 and qi, m (x) = 
((2 — f3) x) e~( 2 ~^ x l 2 (for < (3 < 2, again not identical to previous 
usage of this symbol) . By use of the expansion (|7|) for /3 ^ 1 it follows easily 
that 

rao 

a+1 



qo,n (x) qi, m (x) x a+ dx 



re 

T(a + 2)- 1 / 
J o 

i>oo 

= r (a + 2)" 1 / L« +1 (px) L^ 1 ((2 - (3) x) e~ x x a+1 dx 

(1 _ , r+m ( + jFl (_„, _ m; „ + 2; ^ ' , 

n\ ml \ (I - fly ' 



and the 2F1 -sum equals M n (m; a + 2, (1 - /3) 2 J . Of course this calculation 
is only for illustration and is not novel; the orthogonality of Meixner poly- 
nomials is a simple consequence of that of the Laguerre polynomials and 
the values V (a + 2)~ x / °° (g 0> „ (x)) 2 x a+1 dx = (3- a - 2 ^^ and the similar 
result for qi >m . 

4.1 The transformation matrix 

We consider the analogous computation for the two bases 

(sfl) e- x '^ a+l ^ : m > 0} and {<j> n : n > 0} (from definition |) . The 
scaling for the former is chosen to force equality at m = n = 0. 



Proposition 3 For m,n > let A n = \ y Q+2 1 ra+1 + -^pA , u n = a+ " +1 and 

/ (n, m) = r (a + 2)- 1 f° L« ( X — ) e~ A «* x a+1 dx, 

J \a + 2n + 1 / \a + lj 

then I (0, 0) = (a + l) a+2 ; I (0, n) = = I (n, 0) /or n > 0; A n = (a+1)( 1 1+tt7i) 
and 

Tl \ (I 1 1W1 1 ( a + ^)m-l ( a + / /-i 2\ 

J(n,m) = ((a + l)(l + « n )) ( m -i)!( n -i)i V ~ u n) 

x u™ +n ^ 3 2-F1 (1 - n, 1 - m; a + 3; 1 - u" 2 ) . 

Proof. In / (n, m) change the integration variable to s = A n x and let 

(3 = (A n (a + 2n + l))" 1 = ^gj, then 

I (n, m) = r (a + 2)" 1 A~ a ~ 2 / L° ((3s) L^ +1 ((2 - /?) s) e~ s s a+1 ds. 

Jo 
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By formulae ([/]) and (g) 

a 4 

(n - j) 



3=0 

x ((1-/3) ( a + n + l)-(n-i)) 

= E J ' (a+ f 2 ^. ) , w - 1 - j > " (i - at 1 -' oo • 



for n > while L% (/3s) = 1. By formula (g) 

+ 2 + j) T 
(m-i)! 



m + 2 4- <?'! 

((2 -/3) S ) = E ( 2 " # ^ " ^ L T" ^ 

3=0 



Multiply the two sums and integrate with respect to T (a + 2) 1 e s s a+1 (is 
and by use of the orthogonality relations obtain 

I (n, m) = A~ a ~ 2 (a + 2) m (a + 2) n _ x (-l) m 

x ^ J / /9Qg-2) y 

(„-j)!( m -j)!j!(a + 2) i V (I-/?) 2 7 ' 

for n > 0; and J (0, 0) = = (a + l) a+2 , 1 (0, m) = for m > 0. When 

n > let u n = then A n = ((1 + u n ) (a + l))" 1 ,0 = 1-1^,2-0 = 

l+u n . In the sum replace j by i + 1 (and (n — j)! by (n — 1)!/ (1 — n) i (— l) 4 
and similarly for (m — 1)! ) to obtain the claimed result. ■ 

Certainly the integral could be computed for any scaling of the 
basis, but our particular choice gives the neatest formulae. Motivated by 
the value of I (n, m) we write the following: 

Definition 4 For a > and n = 0, 1, 2, 3 . . . (and u n = a+ "+i ) define the 
function 

h% (x) = u* +1 (x + 1) M n [x; a + 3, u 2 n+1 ) for x > -1. 

Corollary 5 For n > 0, 

I(n,m) - <(„ + 1) (1 + (a + 3) T7' ( °; 2> '- 

ml (n — 1)\ 

x( _i r -i (1 - <)<-*K-X (™ ~ 1) ■ 
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Let 

a+2 (a + 2) m 



e 



-x/(a+l) ^ 



a + 1 



m! 



and recall from Theorem | that ||0„|| 2 = (a + 2n + l) a+3 ^i"" 1 (divid- 
ing by the normalizing T (a + 2)). By the orthogonality of the two bases 
n)_ 

,1/2 



I I (tX, TYl} 1 

the matrix < - — — j- > is orthogonal. The entries in row and 



N ™- ) n,m=0 



column are zero except 7(0,0) = (a + 1) Q+2 = | \<fio\ \ N^ 2 . The row or- 
thogonality demonstrates the orthogonality of the functions {/i"}; indeed 
^Cm^o-^m 1 -^ (n,m) I (l,m) = S n i\ \<p n \ | 2 . After some manipulation this im- 
plies for re, I > 1 that 



E ^r^i (m - 1) fcfL a (m - 1) 
z — ' to! 

m=l 

. (n-1)! ( a + 2 ) 2 (a + n + l) 2a+2 "- ! 

n l- 



(a + 2) n _ 1 ( a + 1) Q+4 (a + 2n + 1) Q+3 n 2 ™" 3 ' 

For n = Z this can be verified from the Meixner polynomial relations: the 
sum equals 

(« + 2 ) E (* + !) "S^n-i (x; a + 3, u\ f ; 

then (x + 1) M n _i (x; a + 3, w 2 ) can be expanded in terms of 
[Mj (x; a + 3, -u 2 ) : j = re — 2, n — 1, n} by use of formula (|iT|), The coeffi- 
cient of M n _i is n<y01 ^^ and together with (||) this produces the desired 
result. The orthogonality of {/i"}^L can also be established by exhibiting 
these functions as eigenfunctions of a self-adjoint operator. 

4.2 A self-adjoint operator 

First let M n _i (x; a + 3; u 2 ) = li^g (x) and substitute in the second-order 
difference equation (JTH) to obtain: 



(x + a + 3) <? (x + 1) + (a + 1) g (x) + xg (x — 1) = (u n + u n x ) (x + 1) # (x) ; 
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replace g (x) by then / satisfies 

X+ x °l + 2 * f + + j^f (*) + / (s - 1) = («„ + n- 1 ) / (*) . 
It is easy to see that the operator 

Vf (x) = X+ x ^ 3 f (x + l)+ - 2) / (*) + / (s - 1) 

defined for functions on x = —1,0,1,2,... such that /(— 1) = is sym- 

00 (a + 2) 1 

metric for the inner product $2) = — -. tt~ /l (x) fi ( x )- Further 

x =o \ x + 1)' 

Vh a = ( U n-±)\ a , f > ^ &nd K ~ ^ = (a + lf ^ 

n_1 u„ n_1 v - " u n n(a + n+l)' 

so here is another proof of the orthogonality of the set {/i" : n > 0}. 
Note that X> can be written as 

2?/(x) = (/(x + l)-2/(x) + /(x-l)) 

+ + -/(*)) + (a+l) (_L- + / (x) , 

x + 2 \x + 1 2 + 2/ 

which is a discrete analogue of the radial Coulomb equation (|l|). It is clear 
that the difference equation 

Vf(x) = -\f(x) (12) 

with initial conditions / (0) = 1 and /(— 1) = has a unique solution for 
any A. To find a formal expression for the solution set 

00 

f(x)=u* (x + l^aji-x^ 
3=0 

with u, a,j to be determined, and substitute in the equation (the series ter- 
minates for each x = 0,1, 2,...). After the usual manipulations we obtain 
(for j > 0): 

-u 2 (j + 1) (j + a + 3) dj+1 
+u (u (a + 2j + 3) + a + 1 - (2 - A) (j + 1)) aj 

- ((u - l) 2 + \u) aj-i = 0. 
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To make this a two-term recurrence, solve u 2 — (2 — A) u + 1 = for u; when 
A < select the solution with < u < 1. Substituting An = — (u — l) 2 in 
the recurrence results in 

-u 2 (j + 1) (j + a + 3) a j+1 + (1 + u) {j (u - 1) + u (a + 2) - 1} a 3 - = 0. 

(13) 

First we handle the special case A = 0, u = 1 with the solution Oj+i = 
(j + ij(j^_a + 3) a j an d denote the corresponding solution / of (12) by where 



7'! 



a hypergeometric series of type i-Fi (and resembling a Laguerre polynomial). 

Second, for u ^ 1 let 7 = M ^ a — - then equation (O) becomes 

u — 1 

(J + 1) (J + a + 3) a j+1 + (u 2 - l) (j - 7) aj = 0, 
with the corresponding solution of ( |T^ ) 

/ ( x ) = u * ( x + l) 2 Fi (-x, -7; a + 3; 1 - u~ 2 ) , 

generally only meaningful for x = — 1, 0, 1, 2, . . . . In this parametrization, 

1 + 7 (a + 1) 2 

and A = — — ; 7-. (The terminating solutions 



a + 2 + 7 ' ( 7 + l)(a + 7 + 2) 

of course were found already, 7 = n for /i°.) The asymptotic behavior of 
Meixner polynomials and their zeros has been investigated by Jin and Wong 
||,||5||, a typical situation being the properties of M n (nx; 7, c) as n — ► 00 for 
fixed x, 7, c. But here c depends on n, and a simple pointwise result is valid: 

Proposition 6 lim h% (x) = (x) for each x € R. 



-2 _ (a+l)(n+2n+3) 

Then ii n +i — > 1 and 



Proof. In h" the argument of the 2-P1 is 1 — n n+1 — . - 2 



-x 



) . (-n) 5 . f ( a + l) ( a + 2 n + 3) V _ + . ^ 



(« + 3), j! V (n + 1) 2 / (« + 3), j! 



for each j as n — > 00. Since 
theorem applies to prove the c 



< 1 the dominated convergence 



aim. 
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For fixed k = 1, 2, 3, . . . the first (in increasing order) k zeros of h% tend 
to the first k zeros of h^. 

In this paper we gave an elementary proof of the strange orthogonality 
of Laguerre functions associated with the wave functions of the Coulomb 
potential and found an analogous orthogonality for Meixner polynomials. 
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